Generalized Lenard Chains, Separation of Variables and Superintegrability 
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INTRODUCTION: SUPERINTEGRABILITY 
AND SEPARATION OF VARIABLES 



The aim of this paper is to establish a new connection 
between the theory of integrable and superintegrable sys- 
tems on one side, and that of bi-Hamiltonian separation 
of variables on the other side. We will provide a the- 
oretical framework for studying separation of variables 
for classical systems, by means of the notion of general- 
ized Lenard (GL) chains. These chains, jointly with a 
couple of compatible Poisson tensors, are the main ge- 
ometrical objects for our bi-Hamiltonian description of 
classical mechanics. These structures guarantee the sep- 
aration of variables in a suitable bi-structured manifold. 

In classical mechanics, superintegrable systems are 
Hamiltonian systems that possess more than N integrals 
of motion functionally independent, globally defined in 
a 2Ar-dimensional phase space (see e.g. [l[ for a mono- 
graph on the topic). These systems are also called non- 
commutatively integrable Q. Especially important 
are the maximally superintegrable ones, i.e. those hav- 
ing 2N — 1 integrals. It turns out that for these sys- 
tems all bounded orbits are closed and the motion is pe- 
riodic The first study in this direction was made by 
Bertrand Q , who derived this result in the case of spher- 
ically symmetric potentials. Among the physically most 
relevant superintegrable potentials are the harmonic os- 
cillator and the Kepler potential, the Calogero-Moser po- 
tential, the Smorodinsky-Winternitz systems, the Euler 
top, etc. [H-[i. 

In quantum mechanics, superintegrable systems are 
also particularly interesting: they possess accidental de- 
generacy of the energy levels. This degeneracy can be 
removed by considering the quantum numbers associated 
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with the additional integrals of motion. Aparadigmatic 
example is offered by the Coulomb atom 0-0 ■ 

Recently, new exarnples of superintegrable systems 
have been discovered In [1], [2Qi], it has been 

proved that the topology of phase space is characterized 
by a bifoliation consisting of an isotropic foliation of in- 
variant tori, and of its coisotropic polar foliation. 

One of the most effective methods to solve Hamiltonian 
systems is to find a complete integral of the corresponding 
Hamilton-Jacobi (HJ) equation through the technique 
of separation of variables. For the sake of clarity, we 
will recall briefly the geometric setting of Hamiltonian 
dynamics (iTJ. 

Let (M, uj) be a symplectic manifold, i.e. a 2n- 
dimensional manifold endowed with a non degenerate 
closed two-form uj, said to be a symplectic form. Such 
a geometrical structure selects a privileged dynamics on 
M, the one given by Hamiltonian vector fields defined by 



= -dH 



{ixn denotes the contraction operator with regard to the 
vector field Xh and d denotes the exterior derivative op- 
erator) or, equivalently 



X 



H 



where uj'' : TM — ^ T* M denotes the fiber bundles iso- 
morphism induced by uj. The function H is said to be 
the Hamiltonian function of the vector field Xh ■ A sym- 
plectic form acting on vector fields is equivalent to a non 
degenerate Poisson bracket defined as 



{F, G] io{XF. Xg) -< dF, Xg > , 



(1) 



(<,> denotes the natural pairing between 1-forms and 
vector fields), i.e. as a skew-symmetric composition law 
on the ring C°°{M) satisfying 

{F,GH} = {F,G}H + {F,H}G 

= {F,{G,H}} + {G,{H,F}}} + {H,{F,G}} 
{F,G} = yF =^ dG ^0. (2) 

A local chart (q,p) ;= (qi, . . . , g„,pi, . . . ,p„) satisfying 
{qi,Pj} = Sij, {qi^Q]} = {Pi^Pj} = is said to be a sys- 
tem of canonical or Darboux coordinates. In such coordi- 
nates the (time-independent) Hamilton-Jacobi equation 
corresponding to a Hamiltonian vector field Xh reads 



dW 



' dqn' 



E 



(3) 



A solution {W, E) with VF(q; oi, . . . , a„), E{ai, . . . , a„) 
and {ai}i<i<n constant parameters, such that 
det[ g^^^ ] 7^ 0, is said to be a complete integral 
of Eq. Q and allows one to solve locally the equation 
of motions for Xh- In fact, a solution of Eq. Q is 



the generating function of a canonical transformation 
that maps the Darboux coordinates (q, p) into a new 
system of Darboux coordinates (q, p) satisfying the 
finite equation of motions 



dE 
oak 



Pk 



(0) 



fc= 1,.. 



(4) 



One can solve the following equation with respect to (q) 



oak 

and by taking into account Eq. one can write down 
the finite equations of motion in the original coordinates. 
In addition, by solving with regard to {ai}i<i<n the fol- 
lowing equations 



dW 



dW , 

Pk = -J, — (q;ai,...,a„ 

oqk 



k — I, . . . ,n (5) 

the function W provides the n involutive integrals pk, 

Pk = flfc = Hk{q, p) fc = 1, . . . , (6) 

whose Hamilton-Jacobi equations admit, by construc- 
tion, the same solution {W,E) than Eq. ([3]). If such 
integrals are globally defined, then Xh is Liouville- 
integrable; we assume this property throughout this 
work. 

Most of the cases in which a complete integral is ex- 
plicitly found, occur when W is an additively separated 
function of the coordinates qi 



W(q; ai, . . . ,a„) 



E 



Wi{qi;ai, ...,an) 



(7) 



In such a case H is said to be separable and the coor- 
dinates (q, p) are said to be separated coordinates with 
regard to H, in order to stress that the possibility of find- 
ing a separated complete integral of Eq. depends on 
the choice of the coordinates. 

One of the classical issues in the theory of separation 
of variables (SoV) is to find criteria to decide if a given 
Hamiltonian function H is separable in an assigned sys- 
tem of canonical coordinates and, in the affirmative case, 
to find a separated complete integral of the Hamilton- 
Jacobi equation. In this regard, a prominent role is 
played by the test by Levi-Civita [22|. It states that 
H is separable in a Darboux chart (q, p) if and only if 
the following n{n — l)/2 conditions. 



= 



dHdH d^H dHdH d^H dH dH d^H 



dqi dqj dpidpj 
dHdH d^H 
dqj dpi dp J dqi ' 



dpi dpj dqidqj 
(1 < i < j < n) 



dqi dpj dpidqj 
(8) 
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are satisfied. 

As we wish to study separable Hamiltonian systems 
that are Liouville-integrable, in principle we can start 
with a set of n independent Hamiltonian functions in 
involution with regard to the Poisson brackets ([T]). In 
this framework, in the tradition of the Italian school, an 
important result has been obtained by Benenti in 23] . It 
gives a characterization of separated coordinates in terms 
of Poisson bracket. 

The Hamiltonian functions {-ffi}i<i<n are separable in 
a set of canonical coordinates {q,p) if and only if they 
are in separable involution, i.e. if and only if they satisfy 



classical Lenard chains, widely known in soliton litera- 
l3C 



31|. 



dH, dHj dH, dHj 



dqk dpk dpk dqk 



, 1 < fc < n, 



(9) 

where no summation over k is understood. However, such 
a theorem as well as the Levi-Civita test are not con- 
structive, since they do not help to find a complete inte- 
gral of the Hamilton- J acobi equation ([3]). In contrast, a 
constructive definition of SoV was given by Sklyanin j24| 
within the framework of Lax systems. 

Definition 1 The Hamiltonian functions {Hi} i<,i<n o.i'^ 
separable in a set of canonical coordinates {q,p) if there 
exist n equations 



^i{qi,Pi;Hi, ...,Hn 







1, 



(10) 



such that det[^^] 7^ 0. They are said to be Sklyanin 

separation equations for {Hi} Ki<n md allow one to con- 
struct a solution (W, E) of the HJ equation In fact, 
solving \1(J\) with regard to pk = , we get 

^ = E J Pk{q'k;H,,...,H,,)\^^^^^dq', (11) 

However, the three above-mentioned criteria of sepa- 
rability are not intrinsic since they require explicit knowl- 
edge of the coordinates (<7,p) in order to be applied. Re- 
cently, a new geometric approach to SoV has been de- 
veloped, based on the bi-Hamiltonian theory ([1^, 26|) 
and on GL chains ([l^-fl^). It has succeeded in giv- 
ing intrinsic and constructive criteria of separability and 
has connected the classical theory of SoV with the mod- 
ern theory by Sklyanin. The bi-Hamiltonian theory of 
SoV is formulated in phase spaces represented by man- 
ifolds endowed with two geometric structures satisfying 
two suitable compatibility conditions. Such structures 
are a symplectic form uj and a Nijcnhuis (or hereditary) 
operator N acting on the tangent bundle of M. For this 
reason such manifolds have been called uiN manifolds. 
Whilst the symplectic form defines the algebra of Hamil- 
tonian vector fields, the Nijenhuis operator defines sets of 
distinguished coordinates that are separated coordinates 
for a special class of Hamiltonian vector fields, those be- 
longing to GL chains, so called as they are extensions of 



ture 

The main result of this paper is the following: we prove 
that, given a generic integrable system on the cotangent 
bundle of the Euclidean plane, the existence of a GL 
chain ensures separation of variables on a ujN manifold. 
We will study explicitly examples of many important 
physical systems, like the Henon-Heiles integrable mod- 
els, the Smorodinsky-Winternitz systems and the Kepler 
potential. For all these systems we will construct explic- 
itly bi-Hamiltonian structures. At the best of our knowl- 
edge, some of the structures we obtain are new. 

Other studies concerning the bi-hamiltonian geometry 
of specific integrable systems that we study in this pa- 
per, are available in the literature [s^-jl^l- A detailed 
comparison between our results and other approaches is 
performed punctually in the subsequent discussions of 
the physical examples. One of the main features of our 
approach is that it provides at once compatible Poisson 
structures and separated coordinates. 

The paper is organized as follows. In Section 2, the 
theory of bi-Hamiltonian manifolds is reviewed. In Sec- 
tion 3, the main geometrical object of our theory, i.e. 
the GL chains, are introduced. A general theorem on 
the separation of variables for systems in T* E2 is pro- 
posed. In Section 4, the previous theory is applied to 
the study of general classes of integrable systems in the 
Euclidean plane. In Section 5, the bi-Hamiltonian geom- 
etry of the integrable Henon-Heiles systems is obtained. 
In Section 6, GL chains are constructed for the classical 
Smorodinsky-Winternitz systems and in particular for 
the Kepler potential. Some conclusions are drawn in the 
final Sec. 7. 



BI-HAMILTONIAN MANIFOLDS AND loN 
MANIFOLDS 

Generally, a Poisson bracket (see, e.g., is defined 
as a skew-symmetric composition law on C°° (M) which 
satisfies only the first two equations of the system ^ and 
not necessarily the third one. Equivalently, it can be de- 
fined by a Poisson bi- vector field, i.e., a skew-symmetric 
linear map P : r*M 1-^ TM such that 

{F,G}p :=< dF,PdG > , 

with a vanishing Schouten bracket 

[P, P]g (a, := Cpp{P)a+P{iPa,dp) a,/3e T*M, 

(12) 

(£ denotes the Lie derivative). In the special case of 
symplectic manifolds, P := (oj^)^^ is a Poisson bi- vector. 
Generalizing ^ the vector field Xq ■— P dG is said to be 
the Hamiltonian vector field with Hamiltonian function 
G. 

Bi-Hamiltonian manifolds were introduced by Magri 
[3^ as models of phase space for soliton equations. 
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Definition 2 A hi-Hamiltonian manifold {M, Pq, Pi) is 

a manifold M endowed with two Poisson bi-vectors fields 
such that 

= 2 [Po, Pi]s (a, /3) Cp„0{Pi)a + Fi(iPoad/3) + 
Cp,p{Po)a + Po{tp,c.df3). (13) 

Such a condition assures that the hnear combination Pi — 
XPq is a Poisson pencil, i.e. it is a Poisson bi-vector 
for each A G C, and therefore the corresponding bracket 
{, — A{, is a pencil of Poisson brackets. Condition 
()13|) is known as the compatibiUty condition between Pq 
and Pi- 

What happens if one of the Poisson tensors, say Pq, is 
invertible, and therefore its inverse is a symplectic opera- 
tor w'' := Pq-^? In this case, the bi-Hamiltonian manifold 
M turns out to be an ujN manifold (see [iOl)- Indeed, 
the composed operator N := PiP^^, thanks to the com- 
patibility condition between Pq and Pi, is a Nijenhuis 
(or hereditary) operator compatible with the symplectic 
form uj. 

Definition 3 A ujN manifold {M,uj,N) is a symplectic 
manifold endowed with an endomorphism of the tangent 
bundle of M , N : TM i— > TM which satisfies the follow- 
ing conditions: 

• its Nijenhuis torsion vanishes identically, i.e. 
\/X, Y e TM 

\NX,NY\ -N{{X,NY\ + {N X ,Y\) + N'^{X ,Y\) = 0; 

(14) 

• it is compatible with u), i.e. the tensor Pi ~ 
N{iJ')^^ is again a Poisson tensor and is compat- 
ible with Pq := (uj^)^^ , according to Definition\^ 



In short, the condition (|T4)) can be rephrased by saying 
that the endomorphism is a Nijenhuis (or hereditary or 
recursion) operator. The adjoint linear map with regard 
to the natural pairing will be denoted by N'^ : T*M ^ 
T*M and will be defined by 

< N^a,X >=< a, NX > . 

The condition iMl) on the endomorphism N, intro- 
duced by Nijenhuis [4l|, has a relevant geometrical mean- 
ing: it implies that the distributions of its eigenvectors 
are integrable according to the Frobenius theorem. Con- 
sequently, under suitable completeness assumption to be 
introduced below, one can select local coordinate charts, 
half of the coordinates being just the eigenvalues of N, 
in which N takes a diagonal form. We suppose that 
at each point x (or in a dense open subset) of M, the 
Nijenhuis tensor field TV admits n distinct eigenvalues 
Xi{x) {i = l....,n) (maximally distinct). Since in a 



generic coN manifold the eigenspaces of N are even- 
dimensional, belonging to the kernel of the skewsimmet- 
ric tensor field Pi — APq, from the above assumption it 
follows that A'' (and the adjoint tensor A'''^) can be put 
in diagonal form. Also, the eigenvalues Xi{x) can be cho- 
sen as coordinate functions in a neighborhood of x, if we 
assume that they are functionally independent in x, i.e. 



dAi(a::)A...AdA„(a;) 7^0 



(15) 



Let a; be a point of an ljN manifold. It will be called 
a regular point if the eigenvalues of N are maximally 
distinct and functionally independent in x. 

From the semicontinuity of the rank function it follows 
that there exist a suitable neighborhood of x whose points 
are regular points. Then we have the following result, 
proved in Let {M,uj,N) be an ojN manifold. 

In a suitable open neighborhood of a regular point, the n 
functions A^ (x) can be completed by quadratures with n 
functions fii{x) such that the chart (A,/x) is a Darboux 
chart for lu and, moreover. 



A^^dA,; = A,:dA, 



A^ d^i = Xid^ii 



(16) 



Then, in a Darboux chart (A, fi) the Nijenhuis tensor A^ 
takes a diagonal form and the coordinates (A, fi) are said 
to be Darboux-Nijenhuis (DN) coordinates. The Poisson 
tensor Pi := A^(a;'')~^ in DN coordinates takes the form 



Pi 



0„ A„ 
-A„ 0„ 



(17) 



where A„ = diag(Ai, • • • , A„) and 0„ is the n x n matrix 
with zero entries. The Darboux-Nijenhuis coordinates 
are just separation coordinates in the bi-Hamiltonian 
theory of SoV. Hereafter, with an abuse of notation, we 
will identify an operator P with its matrix in a suitable 
basis. 

Remark 1 It can be shown that the separated canonical 
transformations 



fi{Xt) , Pi = jj , 



(18) 



with fi a generic invertible smooth function of a single 
coordinate Xi preserve the property ()16p . i.e., 



N^dqi = Xidqi 



N dpi = Xidpi 



(19) 



From a geometrical point of view, we can say that coordi- 
nates (A, /x) and {q,p), related by transformations (|18p 
are adapted to the same coordinate web ]44] - 

GENERALIZED LENARD CHAINS 

After having introduced the geometrical struc- 
tures which define separation coordinates in the bi- 
Hamiltonian theory of SoV, let us characterize the class 
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of Hamiltonian functions which are separable in DN co- 
ordinates. For the sake of concreteness, we wih do this 
in the case of a four-dimensional manifold, warning that 
it has been generalized to a n dimensional manifold .27,] , 

Theorem 4 Let {M,uj,N) be a four-dimensional ujN 
manifold and (Xi, X2, fJ-i, fJ-2) o, DN local chart. Let 
H be a smooth function in M. The DN coordinates 
(Ai, A2, /.ti, /i2) are separated variables for H if and only 
if there exist two smooth functions f and g such that the 
one form 



a = fdH + gN^dH 



(20) 



is an exact one form, i.e., a is the differential of a func- 
tion, say H2 



a = dH2. 



(21) 



In this case, the function H2 is an integral of motion in 
involution with Hi := H and the same DN coordinates 
are separated variables for H2 as well. 

Proof. In the above-mentioned chart, N takes the diag- 
onal form 



N = Ai(-^(g)dAi-l--^(g)d^i)-f A2(T^«)dA2 + T^«)d^2) 
oXi d^i dA2 a^2 

(22) 

Let us suppose that Eqs. (|20p and (pij) are fulfilled. 
Then, it follows that 



dH2 _ fdHi <n\, Ml 



with k = 1,2. Therefore 



{Hi,H2}\k — 



dHi dH2 dHi dH2 
dXk dfik dfik dXk 



^ 



(23) 



(24) 



i.e., Hi and H2 are in separable involution according to 
Benenti's theorem (see formulas (0)), w.r.t. a DN chart. 
Therefore, H2 is an integral of motion for Xh- 

Viceversa, let us suppose that (Ai, A2, /ii, /i2) are sep- 
arated variables for Hi and H2 and let us consider the 
equation 



fdHi + gN^ dHi = dH2 



(25) 



in the unknown functions / and g. In the local chart 
(Ai, A2, ^1, ^2), Eq- (PSj) takes the form 



dH2 

dHi ' 
9^1 



/ + .9A1 



dH2 

/ + 5A2-^, (26) 

dfi2 



We observe that the first equations ([26| and ([27)) coin- 
cide, so do the second equations (p6)) and (|27)) . thanks to 
the conditions (fM)) . Thus the above system of four equa- 
tions reduces to two equations that admit the unique 
solution 



/ - 
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1 



A2-A 



■(A2 



A, 



dH2 

d\i 
2 dHi_ 
d\i 
dH2 



dH2 
\ d\2 \ 

d\2 

Ma 

dX2 
dHi 
d\2 



)■ 



(28) 
(29) 



Then we will say that Hamiltonian functions related 
by Eqs. (f20|) and ([2T|) belong to a GL chain generated by 
{uj,N,H) since, for {f = 0, g = 1), a GL chain reduces 
to a classical Lenard chain. 

Remark 2 We observe that, if (/ ~ — (Ai + A2), g — I), 
we get a quasi-bi-Hamiltonian ( QBH) chain of Pfaffian 
type Wh . \3h l generated by the function H . 

Theorem [4] suggests the following procedure in order 
to classify Hamiltonian systems separable in DN coordi- 
nates: 

1. Choose a Darboux chart (qi, q2,pi,p2) in a 4- 
dimensional symplectic manifold M . 

2. Construct a uN structure which has {qi,q2,pi,P2) 
as DN coordinates. 

3. Search for Hamiltonian function H and for func- 
tions / and g such that make the one form ([20]) 
exact. 

This procedure can be considered as an inverse prob- 
lem, with respect to the direct approach that starts from 
a given Hamiltonian and aims to find separation coordi- 
nates. 

Let us observe that the above method provides the 
integral of motion H2 together with a set of separated 
variables both for H and H2. 



BI-HAMILTONIAN GEOMETRY IN T'^a: 
CONSTRUCTION OF GL CHAINS 

In this section, we wish to apply the procedure previ- 
ously discussed to the study of the bi-Hamiltonian prop- 
erties of systems defined in the cotangent bundle of the 
Euclidean plane. Although we recover well-known re- 
sults, the motivation to study this preliminary case is 
that it allows us to derive the bi-Hamiltonian structures 
for systems separating in one coordinate system in a 
transparent way. The results derived here will be used in 
the following sections to classify multi-separable systems 
and to derive their geometrical properties. Precisely, we 
will study natural Hamiltonian functions 

H — kinetic energy -I- potential energy 
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and we will recover the most general form of the potential 
that makes H separable in Cartesian, polar and parabolic 
coordinates in the Euclidean plane E2- Here M = T*E2. 



Classical Separation of Variables 

Cartesian case 

Let us consider the natural Hamiltonian function 
1 



(30) 



{x , y , Px , Py) being Cartesian coordinates and conjugate 
momenta. According to the requirement (j22p . we choose 
the Nijenhuis tensor Near : T{T*E2) T{T*E2) 

Near = diag{x, y,x,y). 

The one-form ((20|) reads 

a = f {x,y,Px,Py)dH + g{x,y,Px,Py)N^dH, (31) 

and the closure condition da — provides a system of 
nonlinear partial differential equations (PDEs) for /, g 
and V reported in the Appendix, formula (|93|). By com- 
bining Eqs. (j93| . we deduce the interesting differential 
consequence: 



(y - x) gVxy = 0. 



(32) 



The case = is trivial, since it would imply dependence 
of a only on dH. So, if t; ^ and x y, Eq. (15^ implies 
Vxy = 0, i.e. that 

H=l{pl+Pl)+V,{x) + V2{y) (33) 

is the most general Hamiltonian function in T*E2 that 
separates in Cartesian coordinates. In order to get the 
most general form of the integral of motion admitted by 
the Hamiltonian we solve the system ([M)). observing 
that 



k:=?^ + V2{y), 



(34) 



is a particular solution, functionally independent of the 
Hamiltonian (IBSp . Then, we can write down the general 
solution of the system ([M)) as 



/ 



x-y 
F 



F{h,k)+ / Fhdk + G{h) 



x-y 



x^y 



where h := H and F = F {h,k), G = G (h) are arbitrary 
smooth functions of their arguments. The resulting ex- 
pression for the integral of motion is given by a primitive 
function of the one-form (IHTI) 



(35) 



H2 := J F{h,k)dk + J G{h)dh , 

which is independent of Hi := H as 

dHi A dH2 =FdhAdk = Q ^ g = 



In particular, if G = and F = 1 we get the energy 
associated with the coordinate y 



H2^k^^ + V2{y). 



(36) 



Remark 3 Let us observe that g — 1, for any choice of 
f , does not solve Eqs. i93\) . Therefore, neither a Lenard 
chain nor a QBH chain generated by {uj, Near, E[) can 
exist. The same considerations apply to the polar and 
parabolic cases. 



Polar case 



Let us consider the natural Hamiltonian 



H = \U+f2\+y(r^^ 



(37) 



where {r,9,Pr,pe) denote polar coordinates and their 
conjugate momenta. A Nijenhuis operator that separates 
polar coordinates is 



Npoi = diag(r, 6*, r, 6*), 



with 



f{r,0,pr,pe)dH + g (r, 9, pr,pe) N^dH. (38) 

The closure condition for a provides the system (|94l) re- 
ported in the Appendix. By combining these equations 
we get the consequence 



(39) 



gie-r) l^Vr9 + -Vej =0. 
The general solution of Eq. dSH]) { g ^ 0, r ^ 9) is 
V{r,0) = Vi{r) + y^ , 



(40) 



which is the most general potential on E2 that separates 
in polar coordinates. As in the Cartesian case, we note 
that 



k:=§ + V2ie) 



(41) 
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is a particular solution of the system (|94|) . independent 
of the Hamiltonian ([57)) . Then we can state that the 
general solution of the system is 

/ = ^F{h,k)+ J Fhdk + G{h), 



9 = 

t — V 

where h := H and F = F{h,k), G = G (h) are arbi- 
trary smooth functions of their arguments. Therefore, 
the most general expression of the integral of motion 
reads 

H2 := J F{h,k)dk + J G{h)dh, (42) 

which is independent of Hi := H as 

dHi AdH2 ^ Fdh Adk ^0 ^g = 0. 

As a particular case, if G = and = 1 we obtain the 
simple expression of a generalized area integral 



H2 = k = i^ + V2{9). 



Parabolic case 



(43) 



We will study the Hamiltonian 

H=lii4 + V{i,v), (44) 
where ^, 77 are parabolic coordinates given by 

and p^yPri their conjugate momenta. Let us take the 
Nijenhuis tensor 



Npar = diag(^, 77,^,77), 



(45) 



with 



a = fit r],Pi,Pn) dH + g (^, v,p^,Pr,) N^dH. (46) 

The closure condition for a provides the system (^5]) of 
the Appendix. By combining Eqs. (j95p we get the con- 
sequence 

9(V'(){v,.^2!!^)=0. (47) 

The general solution of (|T7)) ( g 7^ 0, ^ 7^ 77) is 



vitv) 



(48) 



which is the most general potential on E2 that separates 
in parabolic coordinates. By following the same proce- 
dure as in the previous two cases, we deduce the general 
solution of system ([M]) 



I = 



F{h, k)+ J Fhdk + G{h) 



F{h,k) ^^r, 



where 



' ^^^^^''+^l($)+V^2(.) 



h 



S.^ + ri^ \ 2 



k ~ 



1 |1P^^^.,,^^^_^.y^^^^ 



with F = F {h, k) and G — G (h) arbitrary smooth func- 
tions of their arguments. Consequently, the most general 
integral of motion admitted by a potential separable in 
parabolic coordinates can be represented as 



H2 := J F{h, k)dk + j G{h)dh (49) 

which is independent on Hi := H as 

dHi A dH2 ^FdhAdk^Q 4=^ g ^ . 
For G = and F = 1 it follows that 

H2 = k. (50) 



INTEGRABLE CUBIC HENON-HEILES 
SYSTEMS 

In this section, we discuss the bi-Hamiltonian geom- 
etry of the integrable cubic Henon-Heiles systems. We 
will show that they are conveniently described by the 
previous theory, and that they admit non trivial bi- 
Hamiltonian structures. The approach we follow is not, 
like the above, an inverse one (i.e. from separation co- 
ordinates toward the Hamiltonian system). Instead, we 
assume the explicit form of the Hamiltonians and con- 
struct the Nijenhuis tensor fields and the corresponding 
chains together with the separation variables. 

For a nice review on integrable Henon-Heiles systems, 
see [4^ and references therein. Besides, new integrable 
perturbations of these systems have been recently ob- 
tained in by means of a Poisson algebra-type ap- 
proach. 

The family of cubic Henon-Heiles systems is defined 
by the Hamiltonian function 



H 



(pI+Py) + 2(^i^^ + '^22;^) + aa;y^-g6a;^ , (51) 
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which is known to be integrable only in three cases: 

(SK) b ^ -a ci^C2^c (52) 
(KdV) b ~ —6a ci , C2 arbitrary (53) 
(KK) 6 = -16a ci = 16c2 (54) 

They correspond, respectively, to stationary reduction of 
the fifth order flow of Sawada-Kotera (SK), Korteweg 
de Vries (KdV), Kaup-Kupershmidt (KK) soliton hier- 
archies [47| . In these tree cases an integral of the motion 
is known 



H. 



(SK) 



PxPy 



1 



a (3 



H. 



(KdV) 



cxy, 



Cl 



(55) 



PyiyPx - XPy) + -(C2 - —)py + C2Xy + 



, (V , 2n 2 C2. Cl 2 

+a(— + ax )y + — C2 - —)y , 
4 a 4 



H. 



(KK) 



Pli^cy^ 



Py 
4 2 2 4 

-a X y 
3 ^ 



Aaxy 
2 



-^ap^Pyy 



ciy 



-a^y^ 
9 ^ 



(56) 



(57) 



The SK cases ([5^ and the KdV case ([5^ have long 
been known to be separable in rotated Cartesian coordi- 
nates [il] and parabolic coordinates [49| , respectively. In 
contrast, only receiitly has the KK case (|54|) been proved 
to be separable in [50|, by means of algebraic-geometric 
methods. 



where 



"-1,1 = (aia4 — a2a3)a; + a2a4(ai — 03)2/, 

ni,2 = a2a4[(ai - 03)2; + (02 - a4)2/)], 

n2,i = -aia3 [(ai - 03)0; + (a2 - 04)1/] , 

n2,2 = aia3(a4 - a2)a; -I- (aia| - a2a3)?/ 



(61) 



The case of the system ([5^ is obtained choosing ai — 
a2 = and 03 = —04 — l/-\/2. The corresponding 

Nijenhuis tensor field is given by 



NsK = 



V2 



X y 

y X 

X y 

y X 



(62) 



the GL chain generated by {ui, Nsk, i?^"^^^) is defined by 

V2 



f = — , 3 = — , y^o, 
y y 



(63) 



and produces the i nteg ral of motion (j55p . 

We recall that in 34^ , a different bi-Hamiltonian struc- 
ture was proposed for the SK case ([52l) . However, being 
constant, it gives no information about separated vari- 
ables. 



KdV— Henon— Heiles system 



SK— Henon— Heiles system 

In order to construct the Nijenhuis tensor field for the 
SK model ([52]) . we shall discuss the general case of Hamil- 
tonian systems separating by means of linear transforma- 
tions of the plane. We have the following result. 

Proposition 1 The most general Nijenhuis tensor field 
associated with systems separating in the normal coordi- 
nates 



Xi^aix + a2y, X2=a3X + a4y, 



aia^ — a2a3 ^ 
(58) 



with associated momenta 
1 



Pxi 

has the form 

N = 



aia^ — 0203 
1 



{a4Px - aspy) , 



{-a2Px + aipy) 



(59) 



0104 — 0203 



^^1,1 »^l,2 

»^2,l »^2,2 

ni^i ^2,1 

ni_2 rt2,2 



(60) 



In the KdV case ([5^ we can put ci = C2 = without 
loss of generality, ci 7^ 0, C2 7^ corresponding to a shift 
of the parabolic web along the y-axis. By transforming 
the Hamiltonian function (|5T|) and ([53|) in parabolic co- 
ordinates we find that it can be described by potential of 
type (gSD with 



V2iv) = -^v' 



The GL chain generated by (lu, N^ar, h[^'^^^ ) provides 

V l<:i=C2=0/ 

the integral ([50]) that, in Cartesian coordinates, coincides 

with dnni). 

Let us observe that the Ni jenh uis tensor field intro- 
duced for the first time in [s^l , |35| in order to construct a 
QBH formulation of the KdV-Henon-Heiles case is noth- 
ing but 



KdV 



(64) 



On the basis of Remark [1] we can state that it defines 
the same parabolic web as (|3S]). However, this result does 
not contradict Remark 3, since the two Nijenhuis tensors 
(Uni) and ([M)) depend on two different realizations of the 
same web. 
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KK— Henon— Heiles system 

In order to construct a Nijenhuis tensor field for the 
KK case (with ci = and a — 1/4), there are (at least) 
two possible procedures. The first one entails the use 
of a well-known canonical transformation between the 
SK and KK cases, that has been introduced and studied 
in [sj, [s^j 51 1- As we have verified, this transformation 



directly maps the Nijenhuis structure of the SK case into 
that of the KK one. 

However, we prefer to follow here a different procedure. 
Precisely, we will consider a mixed problem, in which the 
Hamiltonian of the KK case and its independent integral 
are 



^{KK) 



Pv 



1 2 



r 



Py 
1 

72 



2 2 

Pyxy 



1 3^24 

^PxPyv -j^^y 



y 



(65) 



(66) 



By solving the systems of the 20 differential equations 
[corresponding to conditions i) and ii)] we get the fol- 
lowing structure for the Nijenhuis tensor: 



"1,1 



"-2,1 



"2,2 
"1,4 
"3,2 



fcl 



2 3a;2 jrk^ 

ki^k^x ; h I- — 

ki ki 



<t>{y) 



dX2 



-2x + fca - kok2Py(j)iy) 
koPx<l>{y) , 

u2 



2ki 



Py 

kik2 

W) ' 



I , ,2 3x 2x^3 

1 + fcoKoa: ; 1 ; h k^ 

' Ki ki 



(68) 



where (j){y) is an arbitrary function and kQ,...,ki are 
arbitrary real constants. The four algebraic determining 
equations [condition iii)], which are specific for the chain 
of the KK model, restrict the previous solution and yield 



and the Nijenhuis tensor field is to be determined. In 
our context, the geometrical relevance of the KK system 
is that it does not belong to the Stackel class, i.e. it is 
not separable in orthogonal coordinates in the plane: its 
first integral is a fourth-order one. Moreover, we observe 
that, in contrast to what occurs in all the other examples 
of this paper, the tensor (p7)) is not a complete lift to 
the cotangent bundle (of the configuration space) of a 
torsionless L-endomorphism defined on the configuration 
space 

Therefore, we must assume a general expression for 
the Nijenhuis tensor field in terms of a matrix Nkk G 
A^4x4, instead of a diagonal form. Then, one has to 
impose three conditions: i) that Nkk be compatible with 
the canonical symplectic structure u; ii) that Nkk be 
torsionless; and iii) that Nkk be compatible with the 
GL chain (partially determined by H^^^^ and 'H^2^^^) 
we wish to construct. The algebraic part of the condition 
i) imposes to the Nijenhuis tensor the form 



Nkk = 



"1,1 "1,2 ni^4 

"2,1 "2,2 -"1,4 

"3,2 "1,1 "2,1 

-"3,2 ni^2 "2,2 



(67) 



where n^j are arbitrary functions on T*E2. The remain- 
ing requirements altogether determine a complicated sys- 
tem of 20 nonlinear partial differential equations plus four 
algebraic equations (not reported here). In order to find 
a particular solution of the determining system, a natural 
ansatz is to suppose that the entries of Nkk be rational 
functions of the phase space variables. 



<t>{y) = 

y 

with 



fco = 1, fcl = -, fc2 = 0, ^3=0, ^4=0. 

Then we get the following structure for the GL chain of 
the KK model 



/ 



{xy +&py 



(69) 
(70) 



The eigenvalues of (p7|) 
Ai, 

y" I y" 

together with the conjugate momenta 




Ml,2 



Py3 



-Px + 6^ + 12\± 12 

y y-^ y 



are separated coordinates for the KK-Henon-Heiles sys- 
tem. In fact, they allow to write down the following 
Sklyanin separation equations (jlOp to be written 



pI 



pI 



(KK) 



(71) 
(72) 
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We recall that the above separated coordinates coin- 
cide with those introduced in |5C 



5l| and the Nijenhuis 



tensor field (1571) coincide with a particular case of the 
one introduced in [Hsj, where it has been obtained by a 
completely different method. 



MULTI-SEPARATION OF VARIABLES AND 
SUPERINTEGRABLE SYSTEMS 

In this section, we use the bi-Hamiltonian structures 
constructed in the previous discussion to construct poten- 
tials admitting more than a system of separation coordi- 
nates. In fact, it can achieve that a Hamiltonian function 
belongs to GL chains generated by different and incom- 
patible bi-Hamiltonian structures. In this case, we get 
a Hamiltonian system separable in different coordinate 
system or a multi- separable system together with addi- 
tional integrals of motion that, if they are independent, 
assures superintegrability of the model. Thus, we recover 
in a natural way the Smorodinsky-Winternitz potentials 
in the plane, first discovered in a quantum-mechanical 
context in 54 1, [H^ and studied again in [H^ and 57 1 
from a group theoretical point of view. These are the 



only potentials that are multi-separable in terms of or- 
thogonal coordinates in £2- 

Cartesian and Polar coordinates 

Let us search for the most general potential V{x,y) 
that admits SoV both in Cartesian and in polar coordi- 
nates. To this end, let us write down Eq. p9p (with 
(7 7^ 0) in Cartesian coordinates. It reads 

, } „ {xy{Va;x - Vyy) ~ {x^ ~ y^)Vxy + iyV^ - ixVy) 
V ^ + y 

(73) 

Thus, the potential has to satisfy the system of two PDEs 
V^y = 0, (74) 

xy{V^^ - Vyy) - {x^ - y^)V^y + 3yV^ - 3xVy = 0. (75) 
By substituting the solution of the Eq. (|74l) 

V{x,y) = Vi{x) + V2{y) (76) 
into Eq. ((75|) we get the separated equations 

V(' + = V^' + -Vi = 4a, (77) 
X y 

where a is an arbitrary constant. Their general solution 
is 



n(a;) = -ax + ^, 

V2{y)^\ay^+"-^2- 
2 



(78) 
(79) 



Thus, the general solution of the system (|74p - ([75)) is 

V{x,y) = \a{x^+y^)+'-^ + ^, (80) 

which is nothing but the SWI potential [H^, sum of an 
isotropic elastic potential and an anisotropic Rosochatius 
potential. The Hamiltonian system with the SWI po- 
tential inherits the integral of motion (j36]) from SoV in 
Cartesian coordinates 



H. 



{car)=Py_ a 2 . £2 

2^2^^ 7;2' 



(81) 



and the integral (j43p which, written in Cartesian coordi- 
nates reads 



H. 



(pol) 



i(a:py-2;p,)2 + ci (|)' + C2 (^^^ . (82) 



jicar) 



A simple check shows that the Hamiltonian SWI, H\ 
and H^°^^ are independent. Consequently, the poten- 
tial ([80|) is superintegrable. Finally, we can state that 
the SWI Hamiltonian function generates two GL chains, 
starting from the two incompatible structures (w, Near) 
and {cdjNpoi). Indeed, it can be checked that the two 
Poisson tensors fields 

Pear iV,,,(w^)-i , Ppol NpoiiJ)-^ 

have non vanishing Schouten brackets, i.e., 

[Pear T Ppol] g 7^ ■ 

Cartesian and Parabolic coordinates 



^ 



By repeating the previous strategy, we find that the 
niost general potential V{x,y) that admits SoV both in 
Cartesian and in parabolic coordinates reads 

V{x,y)^aiix^+y^) + cix + ^, (83) 



which is nothing but the SWII potential [55{ . 

The Hamiltonian system with the SWII potential in- 
herits the integral of motion ([36|) from SoV in Cartesian 
coordinates. 



H. 



(car) ^ , a 2 £2 
2 +2'^ 



(84) 



and the integral (HHl) from SoV in parabolic coordinates 
which, written in Cartesian coordinates, reads 

-^2"°'^ = PyiyPcc - xPy) + 2aa:y2 + ^y^ -2c2^. (85) 

^ y 

By checking the functional independence of the Hamil- 
tonian SWII, i/^™''' and H^P°'''\ we obtain again that 
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the potential (|83l) is superintegrable. Moreover, as in the 
previous case, we can state that the SWII Hamiltonian 
function generates two GL chains, starting from the two 
incompatible structures {u),Ncar) and (uj,Npar)- Indeed, 
it can be checked that the Poisson tensors field 

has non vanishing Schouten brackets with Pear 

\Pcar ^ Ppar\ -/- ■ 

In [37I a bi-Hamiltonian formulation for SWI and 
SWII has been proposed. The main difference with re- 
spect to our approach is that, since the two Poisson ten- 
sor fields obtained in [s^l are not compatible [as they do 
not satisfy (fT3|)]. therefore, they do not define a Nijenhuis 
tensor field, and do not provide informations about the 
separated variables admitted by the systems. 

Polar and Parabolic coordinates 

The most general potential y(r, G) that admits SoV 
both in polar and in parabolic coordinates, written down 
in polar coordinates, reads 



r r'- sin 



2Q ' 



(86) 



where 



/3:=J(c3-|) 7:=J(c3 + |-) 

that is the SWIII potential [55]. 

The Hamiltonian system with the SWIII potential in- 
herits the integral of motion (|43| from SoV in polar co- 
ordinates 



(87) 



and the integral (|50l) from SoV in parabolic coordinates, 
which, written down in polar coordinates, reads 



H. 



(par) 



PrSmt 



7 + 2f3cos9 + 7cos^ 



As in the previous cases, the potential (|87p is superinte- 
grable. Furthermore, also the SWIII Hamiltonian func- 
tion generates two GL chains, starting from the two in- 
compatible structures (u],Npoi) and {uj,Npar)- Indeed, it 
can be checked that 

\Ppoli Ppar\g 7^ ■ 



The Kepler model 

The best known example of a Hamiltonian system sep- 
arable both in polar and in parabolic coordinates is the 
Kepler model 



H 



1 



K 



{pI+pD 



Here, we recover from the theory illustrated above a new 
bi-Hamiltonian structure and a GL chain for this model. 
Precisely, by exploiting the freedom in the choice of sep- 
arated coordinates discussed in Remark [U we take as 
a Nijenhuis tensor field that separates polar coordinates 
the following one, 

NpoiK — diag(r, tan 6*, r, tan 6*), 

that in Cartesian coordinates reads 



N, 



polK 



ni,i 










ni,2 













"-3,2 




"1,2 


-n3,2 





"1,2 


"2,2 



(89) 



where 



"1,1 

"1,2 
"2,2 
"3.2 



X 

r 



_T_ xy_ 

j,2 J, '' 



xy_ 



r 



I y y. 
(— - -)Px 



-_y_ 

r2 



X 

-)Pv ' 



(90) 



with r 



{uj,Npoik,Hk) with 
y — XI 



The GL chain generated by 



.9 = 



y 



provides as a second integral of motion the square of the 
modulus of the angular momentum: 

Furthermore, let us take as a Nijenhuis tensor field that 
generates the parabolic web the following tensor. 



NparK = diag(^'^,-Tf,i'^,-rj^) 



(91) 



which in Cartesian coordinates has the linear represen- 
tation 



parK 



2x y 

y 

py 2x y 

-Py y 



(92) 
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The GL chain generated by (w, NparK, Hx) with 

/ = 2a;, g = -l, 

provides as a second integral of motion the x-component 
of the Laplace-Runge-Lenz vector 



H. 



{parK) 



Py{xpy - yp^ 



Surprisingly, the Kepler Nijenhuis tensor (1911) coincides 
with the Nijenhuis tensor (j64p . constructed by other au- 
thors for the KdV-Henon-Heiles system. This is due to 
the fact that NparK generates the same parabolic web in 
which both systems are separable. 

In [3^ a different bi-Hamiltonian formulation for the 
three-dimensional Kepler problem has been introduced. 
The main motivation of the authors was to prove that 
the existence of a recursion operator does not necessar- 
ily provide additional conservation laws. The structures 
obtained in 33] are expressed in terms of action-angle 
variables, and give no information about the separated 
variables. 

In ^3|] another recursion operator for the Kepler po- 
tential has been proposed. However, the difference with 
respect to our approach is that their recursion operator is 
not compatible with the canonical symplectic structure. 



FUTURE PERSPECTIVES 

In this work, we have proposed a general formalism for 
treating the geometry of both intcgrable and superinte- 
grable systems on a bi-Hamiltonian setting. The present 
approach for the sake of concreteness has been formu- 
lated in the Euclidean plane. However, there is no the- 
oretical restriction in extending it to higher-dimensional 
cases. Also, it seems interesting to include in the present 
analysis integrable and superintegrable systems defined 
in curved spaces. It would be very interesting to derive 
a quantum formulation of the present theory. 



Appendix 

Here we report the explicit expressions of the systems 
of differential equations quoted in Sec. HI. 
a) Cartesian system. 

Pyifp. + ygpJ-Pxifpy + xgpj = 0, 
Pxifx + xg^ + g) - {fp^ + xgpjVx = 0, 
Pyifx + ygx) - ifpy + xgpjVx = 0, 

Pxify + xQy) - {fp^ + ygpJVy = 0, 

Pyifv + ygy + g) - (/p„ + 2/5pJK/ = 0, 

{x - y)gVxy + ify + xgy)V^ - {f^ + ygx)Vy = 0;.93) 



b) Polar system. 
Po 

^ (/p. + %J - Pr {fpe + rgp,) = 0, 
Pr ife + rgg) - {fp^ + 9gp^) Vg = 0, 
5 {fe + ege + g) - (Jp, + egp,)Vg = 0, 
P^ 

Pr ifr + rgr + g) + {fp,. + rgp,.) - ifp^ + rgp^) Vr = 0, 
5 (/. + Ogr) + ^2g {r-e) + ^ {fp, + rgp, ) 

- ifpe +rgpe) Vr = 0, 

{6 - r) gVrg + [fr + Vg - {fe + rgg) Vr 

+^(/e + r.ge) = 0. (94) 

c) Parabolic system. 

Pi {fp„ + ^gp„) ~ Pv {fps + Wpe) = 0, 

Pi — 7^-rz72 — + I ^77^— ;t2 - ^« (/p£ + ^5pJ = 0, 



U + ^g) 



Pv 



X (/p„+e5p„) =0, 



f + vg 



2?7P{ 



(/ + vg) 



Pi 



f + ^g 

^2 _^^2 



X (/p^ +vgp() = 0, 
ifv + vgy+g) , f pI+Pv 



(^2+^2) 

g{^~v) Vir, + if, + ^gn) Vi - {fi + m) Vr, + {pI + p^) 



vlif + vg) 



Pj+Pr, 

ie+v')\ 



0. 



(95) 
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